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Abstract 

We present results concerning resolvent estimates for the linear op- 
erator associated with the system of differential equations governing 2 
dimensional perturbations of plane Couette flow. We prove estimates on 
the 1/2 norm of the resolvent of this operator showing this norm to be 
proportional to the Reynolds number R for a region of the unstable half 
plane. For the remaining region, we show that the problem can be re- 
duced to estimating the solution of a homogeneous ordinary differential 
equation with non-homogeneous boundary conditions. Numerical approx- 
imations indicate that norm of the resolvent is proportional to R in the 
whole region of interest. 



1 Introduction 



We consider the Initial Boundary Value problem 
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H 



(1.1) 



w t + (w- V)W + [W ■ V)w + Vp 
V • w = 

w(x,0,t) = w(x, l,t) = (0,0) 
w(x, y, t) = w(x + 1, y, t) 
w(x,y,0) = (0,0) 

where w : Rx [0,1] x [0, oo) — > M 2 is the unknown function w(x,y,t) — 
(u(x, y, t), v(x, y,t)). W is The vector field W = (y, 0) and the Reynolds number 
R is a positive parameter. The forcing H(x,y,t) = (F(x,y,t),G(x,y,t)) is a 

/>oo 

given C°° function satisfying satisfying / ■, t)\\ 2 dt < oo and V • H = 

Jo 

for all (x,y,t) G R x [0, 1] x [0, oo). These equations are the linearization of 
the equations governing 2 dimensional perturbations w(x,y,t) of W = (y,0), 
known as Couette flow, which is a steady solution of 

U t + (U- V)E7 + VP = —AU 
R 

V • U = 

U(x,0,t) = (0,0) ( L2 ) 
U(x,l,t) = (1,0) 
U(x,y,t) = U(x+l,y,t). 
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The pressure term p(x,y,t) in (1.1) is determined in terms of w by the linear 
elliptic equation 



Ap = -V ■ ((w ■ V)W) - V • ((W ■ V)w) 
Py(x,0,t) = -^v yy (x,0,t) 

Py(X,l,t) = -jjVyy(X,l,t). 



(1.3) 



We note that p depends linearly on w, and is determined up to a constant. The 
estimates derived in this paper are independent of p. With p given by the above 
equation, the solution w of (1.1) remains divergence free for all t > 0. Therefore 
we drop the continuity equation and write the problem as 

f w t = ^Aw - (w ■ X7)W - (W ■ V)w -X7p + H =: Crw +H ^ ^ 
1 w(x,y,0) = (0,0) 



with w satisfying the boundary conditions described in (1.1). Note that Cr is 
a linear operator on ^(ri), f2 = [0, 1] x [0, 1]. The L2 inner product and norm 
over O are denoted by 



(wi,w 2 ) = / wi ■ w 2 dxdy ; || w|| 2 = (w, w). 



As motivation, consider a general linear evolution equation 



w t — Cw + H 
w(x,y,0) = 0. 



(1.5) 



Formally, after Laplace transformation with respect to t, the transformed equa- 
tion is sw = Cw + H for Re (s) > 0, where w is the Laplace transform of 10. 
The solution of the transformed equation is given by 



= {sI-L)- 1 H. 



(1.6) 



wherel is the identity operator, and then we get the solution of (1.5) by applying 
the inverse Laplace transform to w. Moreover, from (1 



|^ s )|| 2 <||( s z-£)-Tl|ffM 



(1.7) 



We recall the following definition: 

Definition 1 Let C be a linear operator in a Banach space X and X the identity 
operator. The resolvent set p(C) C C of C is the set of all complex numbers s 
such that the operator (sX — C)^ 1 exists, is bounded and has a dense domain in 
X . The linear operator (sT — C)^ 1 , for s 6 p(£), is said to be the resolvent of 
C. The set o~(£) = C \ p(C) is the spectrum of C. 
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Therefore, if the complex number s is such that s G p(C), the formal argu- 
ment above to express the transformed function w in terms of the transformed 
forcing H by (1.6) is valid and we have the estimate ([O]). According to Ro- 



manov |Q, the unstable half plane Res > is contained in the resolvent set 
of the linear operator Cr defined in ( |l.4| ) , for all positive Reynolds numbers 
R, and the eigenvalue of Cr with largest real part is at a distance from the 

imaginary axis proportional to — . Our aim in this paper is to estimate the 

R 

resolvent constant sup \\{sl-C R )- x \\. For each R, we derive estimates on the 

Re s>0 

1/2 norm of the resolvent for the region \s\ > 2\/2(l + a/R)- For the remaining 
region, we prove that we can reduce the problem to estimating the norm of the 
solution of a homogeneous ordinary differential equation with non-homogeneous 
boundary conditions. This is the main contribution of the paper, since in prin- 
ciple it simplifies the problem either if one wants to prove the estimates ana- 
lytically or use numerical computations. We perform numerical computations 
that indicate that the Li norm of the resolvent is proportional to R also for 
\s\ < 2y / 2(l + \/R)- We note that, according to Romanov this is the best 
possible dependence of the resolvent on the Reynolds number. This dependence 
is better than the 3 dimensional case, since results in Liefvendahl & Kreiss||, 
Reddy & HenningsonQ, Trefethen et aZ. ||, Kreiss et aZ.Q, also partly based 
in computations, indicate the L 2 norm of the resolvent to be proportional to 
R 2 . In all the papers above, the computations are performed using different 
methods than the one used here. 



2 Resolvent estimates 
2.1 s bounded away from 

Theorem 1 If\s\ > 2\/2 f 1 + Vr) , then \\ (si- Cr.)^^ < -^(l + Vfi) 2 < 1. 



Proof: After Laplace transformation in t, the first equation in (1.1) is trans- 
formed to 



sw = — Aw 
R 



(w ■ W)W - (W ■ V)w -Vp + H = C R w + H 



(2.1) 



with w satisfying the boundary co ndit ions in the space variables described in 
( |l.l[ ). Taking the inner product of (2.1) with w, we have 

(w,sw) = (w, ^-Aw)-(w, (w-W)W)~(w, (W-W)w)-(w,Wp) + (w,H). (2.2) 

R 

Integrating by parts and using the divergence free and boundary conditions, one 
can prove that (w, Vp) = and (w, (W ■ \7)w) is purely imaginary. Therefore, 



S IH| 2 + oll-D^II 2 = -(w,(w ■ V)W) - (w,(W-V)w) + (w,H) =: P (2.3) 
it 
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where Dm denotes the derivative of w with respect to the space variables, and 
its Frobenius norm. P satisfies 

\P\ < \\w\\\\(w ■ V)W\\ + \\w\\\\(W ■ V)w\\ +IHIHFH (2.4) 

which implies, since ||£W|| = 1, 

|P|<H| 2 + \H\\\Dw\\ + |H|||#||. (2.5) 

Since (w, (W ■ V)wJ) is purely imaginary, taking the real part of ( |2.3| ) gives 

Res||w|| 2 + — ]|ZWo|| 2 < \(w,(w- V)W)| + \{w,S)\ < \\w\\ 2 + \\w\\\\H\\ (2.6) 



R 



and therefore 



Res||w|| 2 < ||w|| 2 + 11^1111^11 



(2.7) 



Note that (2.6) and fl2.7|) are valid for all s satisfying Res > 0. To complete the 
proof, we consider two separate cases: 

Case 1: Res > |Ims|. 



In this case, Res > — p, and we choose s such that — -= > 1. Using (2.7), 
V 2 2v 2 



2V^2 

which implies 



V2 2V2 

N 2 N| 2 <8||#|| 2 



||w|| 2 < (Res - l)||w|| 2 < 



(2.8) 



Case 2: |Ims| > Res > 0. 



In this case, |Ims| > Taking the imaginary part of equation (2.3) gives 
v2 



By <M), 



and then 



|Ims||H| 2 < \P\ < \\w\\ 2 + \\w\\\\Dw\\ + \\w\\\\H\\ 



|||^|| 2 <||^|| 2 + ||^||||ff||<(||^|| + ||ff|r 2 



\\Dw\\ < VR,(\\w\\ + \\H\ 



Using this estimate in (2.£), we get 

|Ims|||w|| 2 < (1 + VR)\\w\\ 2 + (1 + V^K)|]w|||j-H"|| 
Isl 

Choosing s such that — j= > 1 + \/~R, we have 



2V2 



2V2 1 



>\\ 2 < (|Ims| - (l + VR)) \\w\\ 2 < (1 + VR)\\w\\\\H\ 



(2.9) 



(2.10) 
(2.11) 
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and this implies 

|s| 2 ||w|| 2 <8(l + Vi?) 2 ||i/|| 2 . 
From the two cases, we conclude that 

|s| > 2V2(1 + VR) =► |s| 2 ||w|| 
which implies the desired estimate 

\\( s I-£ R )-if<*(l + jRf<l, 



<8(1 + VR) 2 \\H\\ 



(2.12) 
(2.13) 

(2.14) 



valid in the region |s| > 2-\/2(l + V~R) of the unstable half plane Res > 0. 



2.2 s<2V2(l + VR 



For this case, write the problem (1.1) componentwise: 

1 



-Au + F 



R 

Av + G 



ut + yu x + v + p x 

] 

vt + yv x +Py = -j 

U x + Vy = 

u(x, 0, t) — u(x, 1, t) — v(x, 0, t) ~ v(x, 1, t) = 
u(x,y,t) = u(x + l,y,t) 
v(x,y,t) = v(x + l,y,t) 
I u(x,y,0) = v(x,y,0) = 



(2.15) 



Introduce the stream function ip by 



tPy = -U. 



(2.16) 



Laplace transform problem ( 2.15 ) with respect to t, expand in a Fourier series in 
the periodic direction x. The equations for the transformed functions u(k, y, s), 
v(k, y, s), p(k, y, s), F(k, y, s), G(k, y, s) are 



su + ikyu + v + ?fcp = u + — u 

R " 
1 



+ ikyv + p y = -—v + —v yy 
R 



iku + v v 



0. 



+ G 



F 



(2.17) 



Differentiating the first equation in (2.15) with respect to y, the second with 
respect to x and subtracting the second from the first, we eliminate the pressure 
p, and using the divergence free condition, we see that the transformed stream 
function ip satisfies the following boundary value problem for an ordinary dif- 
ferential equation with three parameters s, R, k: 



R 



2 4 + ^y)^" + ( sk2 + ^ + ^y)^ 1 (2.i8) 



ip{k, 0, s) = ip{k, 1, s) = V'(fc, 0, s) = i>'(k, 1, s) = 
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where I :— F y — ikG and ' denotes the derivative with respect to y. We also 
d 

use the notation D = — , and write the problem as 
oy 

TT oV 7 = (i^ - (a + | + iky)) (P 2 - = I ^ 
i/>(k, 0, s) = i/j(k, 1, s) = i/j y (k,0, s) = i) y (k, 1, s) = 

where T = —T> — [ s + — + iky and Tq =T> — k are differential operators 
R \ R J 

depending on the parameters R, k and s. To derive the resolvent estimates, we 

use the following Lemma: 

Lemma 2 If for all k S Z and /or all s <E C such that \s\ < 2\[2{\ + i/ie 
solution ip(k,y,s) of (2. IS) satisfies 



k 2 U(k,; S )\\ 2 <CR 2 (\\F(k,; S )\\ 2 + \\G(k 1 - 1 s)\\ 2 ) 

U'(k,-,s)\\ 2 <CR 2 (\\F(k,-,s)f + \\G(k,-,s)f) 



(2.20) 



where C is a constant independent of s, R, k, F , G, then 

\\(sl - Cr^W 2 < 2CR 2 

for \s\ < 2y/2(l + y/R). 

Proof: If for all k G Z and for all s £ C such that |s| < 2^(1 + 

I U'(k,-,s)f<CR 2 (\\F(k,-,s)\\ 2 + \\G(k,-,s)\\ 2 ) { - ' 



then by the definition (2.16) of the stream function, 

\ ||«(fc,., s )|| 2 = ||^(fc,., s )|| 2 <Ci? 2 (||%,., s )|| 2 + ||G(fc,-, s )|| 2 ) 
and since 

oo 

iiu(.,.,.)n 2 = J2 ii^-' s )ii 2 

k— — oo 

oo 

||5I(,, S )|| 2 = £ \\v(k,;S)\\ 2 
k—-—oo 

||%-, S )|| 2 + ||G(-,-, S )|| 2 = J2 (||F(fc,-, S )|| 2 + ||G(fc,-, S )|| 2 ), 



(2.22) 



( p.22j ) implies 

/ P(-,-, s )|| 2 <G J R 2 (||i?(.,., s )|| 2 + ||G(-,-, s )|| 2 ) 
|R,, S )|| 2 <Gi? 2 (||F(,, s )|| 2 + ||G(,, S )|| 2 ) 
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and this implies || (si - Cr)' 1 f < 2CR 2 for \s\ < 2^2(1 + y/R). 

This Lemma shows that to estimate the norm of the resolvent, it is enough 
to prove estimates of the form (2.21) for the solution ip of Q2.18 ). To prove those 
estimates for ip, take the inner product of the differential equation in (2.18) with 
tp and obtain 

i ok 2 - A- 4 - 

-(^^"")-(^(s+ — +iky)^") + (^(sk 2 +-+ik 3 y)^) = (V.,/). (2.24) 



Through integration by parts, 



^,i,"") = \W'f 



(2.25) 



Therefore equation ( [2.24 ) becomes 



^H"\\ 2 +( + s ) W\\ 2 +[ sk 2 +'—) yj\\ 2 -ik(TP,yiP")+ik 3 (iP,y4>) = 



7? 



R 



R 



Again through integration by parts, we have (ip,ytp") = — (ip,ip') — (ytp',ip'}, 
and since ip is a real function, (ip,ip') is purely imaginary, (yip',tp') and (ip,yip) 
are real numbers. Therefore, taking the real part of the previous equation we 
get 



1 

R 



2P 
R 



+ Res) \\tp'\\ 2 +[ (Res)fc 2 + ^- ) = Re{ip,I)-Re(ik{tp,i)')) 



R 



and since Re s > , this equation implies 
1 



R 



'II 1 1 2 



-^U'f + ^Uf-\H\(i'^')\<\(^i)\ 



(2.26) 




We separate the analysis into three different cases for k € Z: |fc| > w — -= 
\ R 

< |fe| < y — ^= and the special case fc = 0. 



Case 1: |fc| > 




For this case, we prove the following: 



R 

Theorem 3 If \k\ > then 
\J2 



k 2 \$(k, ,s)|| 2 < l6R 2 (\\F(k, s)|| 2 + \\G(k, s)\\ 2 ) 
U'(k, •, S )|| 2 < 16i? 2 (||F(fc, s )|| 2 + \\G(k, ;.s)\\ 2 ). 
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Proof: Since \(ip,^')\ < 
implies 



< 



— I 



I 2 + ^ll?H 2 , inequality (p|) 



2fc 2 



R R 



R 



(2.27) 



and since |fc| > 



i? k 4 R k 4 

— = , we have — > — - and (2.27) implies 

V2 R 4 2R vc — " 



^ll?'ll 2 + ^ll?ll 2 + ^ll^ll 2 <l(V",/)|. 



(2.28) 



Since equation ( 2.18 ) is linear and the forcing is I — F y — ikG, it suffices, to 
get the desired estimates for the solution ip } to prove estimates for tpi and ip2, 
solutions of the boundary value problems 



1 



R 



2k 2 



^»- {s + - + iky)^ +[ se + - + ik*y)^ = F y 



k 4 

.2 , ft , „-i.3„ 



R 



(2.29) 



V>i [k, 0, s) = 1, «) = 0, s) = $ (A, 1, s) = 



and 



1 



R 



2k 2 



R 



-V>2 - s + — + iky ^2 1 + sfc + s + ik V V-2 = ikG 



R 



(2.30) 



V- 2 (fc, 0, s) = ^ 2 (fc, 1, s) = V4(fc, 0, a) = ip' 2 {k, 1, s) - 0, 
since V> is given by ip = ipi + V>2- 
Estimates for ipi- Through integration by parts, 

\(M)\ = \$ 1 ,F)\. 



Using the Cauchy-Schwarz inequality, ( J2.28| ) with forcing F y yields 
^ll^l| 2 + ^ll^l| 2 + |^ll^l| 2 <ll^llll^ll 



(2.31) 



and since all the terms on the left hand side of this equation are positive, we 
have 

fc2 -WH 2 < 11^1111^11 => k 4 U[\\ 2 < R 2 \\F\\ 2 (2.32) 



and 



and 



R 
2R 

iii^ir<ii^iiii^iK* 2 iiw<^ii 2 - 



||^i|| 2 < Mllll^l ^ 6 ||Vi|| 2 <2i? 2 ||i?|| 



(2.33) 
(2.34) 
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Estimates for ip2' By ( 2.28) ) with forcing ikG and the Cauchy-Schwarz in- 
equality, we have 

^ll^f + j^W + ^\\M 2 < \k\U 2 \\\\G\\ (2.35) 

and then, similarly to the argument used to prove the estimates for %p\ above, 
we conclude 

' fc 6 ||^ 2 || 2 <4i? 2 ||G|| 2 
fc 4 ||^|| 2 <2i? 2 ||G|| 2 (2.36) 
fc 2 |fe'|| 2 < 2i? 2 ||G|| 2 . 



Therefore tj), the solution of problem (2.18), satisfies 



k 2 U"\\ 2 + k 4 \\ip'\\ 2 + k e \\ip\\ 2 < 16i? 2 (||F|| 2 + ||G|| 2 ) = 16i? 2 ||#|| 2 . 
Case 2: < Iftl < 

V ^ 

For this case, we prove that we can reduce the problem to the study of the 
solutions of a homogenous 4th order ordinary differential equation with non- 
homogenous boundary conditions. We compute the norms of those functions 
numerically. The computations indicate that the norm of the solutions of these 
simplified problems do not grow as R grows. This implies the norm of the 
resolvent of the operator Cr to be proportional to R. We begin by noting that 
we can restrict ourselves to the case when s = i£ is purely imaginary. This is 
a consequence of the following theorem for holomorphic mappings in Banach 
spaces (Chae [0): 

Theorem 4 (Maximum Modulus Theorem) Let U be a connected open subset 
of C and f : U — ► E a holomorphic mapping, where E is a Banach space. If 
\\f(z)\\ has a maximum at a point in U , then ||/(z)|| is constant on U . 

Applying the theorem to the holomorphic function f{s) = (sX — Cr,)^ 1 defined 
on Re s > 0, taking values in the Banach space of bounded linear operators on 



Li2{£l) and noting that ( 2.14 ) implies lim ||/(s)|| = 0, we conclude 

|s| — >oo 

sup H&Z"— £r|| = sup||i£Z"— Cr\\. (2.37) 

_Res>0 CeK 

We note that since we just need to consider s — i£, £ S R, and in this case 
l s l = l£l = |Ini s|, the proof of Theorem |l| is valid for \s\ — \lms\ > 2(1 + Vi?). 
This restricts a little bit the parameter range for the numerical calculations. 
We now reduce our problem to the study of solutions of a homogenous ordinary 
differential equation. To this end, first consider the second order system 

/ 1 1,2 \ _ 

Th = ( —V 2 - (s+ — + iky)j h = I = F y - ikG , h(0) = h(l) = 
T a g - (V 2 - k 2 )g = h , 5 (0) = 5 (1) = 
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for s = i£, £ G R. Taking the inner product with h in the first equation and with 
g in the second, and since the boundary conditions for both equations imply 
that the boundary terms after integration by parts vanish, we get 



^"i 2 + fc 2 ||/ l || 2 <G 2 i? 2 (||F|| 2 + ||G|| 2 ) 



(2.38) 



where Gi, Ci are constants independent of i?, s, k, F, G. Combining those two 
inequalities, it follows that 



k 2 \\g"\\ 2 + k 4 \\g'\\ 2 + k 6 \\g\\ 2 < Cii 2 (||F|| 2 + ||G|| 2 ) 



(2.39) 



where G is a constant independent of R, s, k, F, G. Note that g satisfies 



TT Q g = I — O - 
g(0) = g(l) = 0. 



- + iky ) ) (V 2 - k z )g = F y - ikG 



(2.40) 



Let g'(0) = a and </(l) = /3. The numbers a and (3 can be estimated using a 
1-dimensional Sobo lev type inequality. Indeed, Ig'l^ < || <?'|| 2 + \\g'\\ \\g"\\ . Using 
the estimates ( 2.39| ) and that fceZ satisfies 1 < we conclude that 



| 5 U<Ci? 2 (||F|| 2 + ||G|| 



(2.41) 



and therefore 



\a\ 2 = \g'(0)\ 2 < Gi? 2 (||F|| 2 + ||G|| 2 ) 
\f3\ 2 = \g'(l)\ 2 <CR 2 (\\F\\ 2 + \\G\\ 2 ), 



(2.42) 



representing again by G a constant independent of s, R, k, F, G. Now, let 
S(k, y, s) be the solution of the problem 



1 



R 

(5(0) = 6(1) = 
8'(0) = a 
{ 5>(l)=p. 



k 



TT 5 = -V 2 - s+ — + iky) (V 2 - k 2 )5 = 



R 



(2.43) 



Then, ip(k, y, s) — g(k, y, s) — S(k, y, s) is the solution of ( 2.19|) . Indeed, 
TT ^ = TT {g -6)= TT g - TT 5 = F y — ikG 



and 



$(k, 0, s) = g(k, 0, s) - 5(k, 0, s) = 
$(k,l,s) =g(k,l,s)-6(k,l,s) = 
$'(k, 0, s) = g'(k, 0, s) - 5'(k, 0, s) = a - a = 
f(k, 1, s) = g'(k, 1, s) - S'(k, 1, s) = $ - = 0. 
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Since we already have suitable estimates for g , our problem is reduced to 
deriving estimates for the solution of 



TTqS = (j^D 2 - 

5(0) = 5(1) = 
5'(0) = a 



where |of < CR 2 (\\F\\ 



R 



— - + iky)) (V 2 - k 2 )5 = 



(2.44) 



||G|| 2 ) and \(3\ 2 < CR 2 ( 
it remains to estimate k 2 \\5(k 1 ■, s)\\ 2 and ||<5'(fe, •, s)| 



!|F|| 2 + ||G|| 2 ). Therefore, 
2 for the parameter range 



€ Z X E, 1 < \k\ < . —= , |^| < 2(1 + y/R), where s = i£. To study this 
V v2 

problem numerically, we simplify it as follows: if 5% is the solution of 



and 52 is the solution of 



TT 6i = 
d 1 (0)=5 1 (l) 
5[(0) = 1 
51(1) =0 



TT 6 2 = 
5 2 (0) = 5 2 (1) 
5' 2 (0) = 
= 1. 



= 



(2.45) 



= 



(2.46) 



then 5 = a8\ + (35 2 is the solution of ( |2.44| ) • Therefore 

fc 2 ||<5(fc, •, s )|| 2 < 2|a| 2 fc 2 || ( 5 1 (fc, •, s )|| 2 + 2|/3| 2 fc 2 ||<5 2 (fc, •, s )|| 2 

and 

||<5'(fc,-, S )|| 2 <2|a| 2 ||^(fc,-,s)|| 2 + 2|/?| 2 ||^(fc,., S )|| 2 . 

Thus we can restrict ourselves to the dependence on R of k 2 \\5j(k, ■, s)\\ 2 and 
||5j(fc, •, s)|| 2 , for j — 1, 2. Moreover, since \\ip(k, •, s)|| 2 = ||^(— k, •, — s)|| 2 , where 
ijj is the solution of (2.18), we can restrict ourselves to < £ < 2(1 + VR). We 
performed numerical computations using the MATLAB 6.0 built-in boundary 

/ R 

value problem solver BVP4C, for the parameter range 1 < < a —=, < 

V v 2 

C < 2(1 + and values of R from 1 up to 10000. For £, we used a mesh 
with variable number of points for each R. BVP4C makes use of a collocation 
method, and we performed computations with different absolute and relative 
tolerances, using continuation in the Reynolds number for the initial guess of 
the solution. The results were similar for all cases. Therefore, even though the 
problem is stiff for some of the parameter values, the results shown in figures (|l|), 
(||), (||), ([IJ) in the following pages should be reliable. For different values of the 



12 



Resolvent estimates for 2D perturbations of plane Couette flow 



Reynolds number R, we plot the maximum of k 2 \\Sj(k, •, s)\\ 2 and \\S' (k, •, s)|| 2 
/ R 

for 1 < k < \i—j=, < £ < 2(1 + VR). The results indicate that 
V v2 



fc 2 ||(5,(fc,-, S )|| 2 <l and ||£'(fc,., S )|| 2 <l i = 1,2. 



(2.47) 



Therefore, 

fc 2 ||5(fc,, S )|| 2 <2H 2 fc 2 ||^fc,, S )|^ 
that is, 



and similarly 



k 2 \\S(k,; S )\\ 2 <CR 2 (\\F(k,; S )\\ 2 + \\G(k,; S ) 



\\5'(k,;s)\\ 2 <CR 2 (\\F(k,;s)\\ 2 + \\G(k,; S )\\ 2 ) 



(2.48) 



(2.49) 



/ R ~ 
for 1 < |fc| < 4 / — p , |C| < 2(1 + VR), where G is a constant independent of R, 
V v2 

s, fc, F, G. Those inequalities imply 

f fc 2 ||^(fc,-, S )|| 2 <C J R 2 (||F(fc,-, S )|| 2 + ||G(fc,-, S )|| 2 ) 
\ ||?(fc, S )|| 2 < CR 2 (\\F(k, s)|| 2 + ||G(fc, S )|| 2 ), 

where G is a constant independent of R, s, k, F, G. Those are the desired 
estimates for ip. 



Figure 1: max/c 2 ||(5i(fc, •, ^)l| 2 for 1 < A: < -2(1 + VR) < £, < 2(1 + s/R). 
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Figure 2: max UUk, i£)f for 1 < k < ./X, -2(1 + < £ < 2(1 + \/#). 

fc,£ V 
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Figure 3: maxk 2 \\5 2 (k, -,iOf for 1 < k < x f-%, -2(1 + VI) < £ < 2(1 + VK). 
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Resolvent estimates for 2D perturbations of plane Couette flow 



Figure 4: max \\S' 2 (k, ■, i£)\\ 2 for 1 <k< J^=, -2(1 + ^R) <£ < 2(1 + 77?). 
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Case 3: k = 

For this case, the differential equation for ip is reduced to 



(2.50) 



^(0, 0, a) = ^(0, 1, s) = V'(0, 0, a) = ^(0, 1, s) = 0. 
By an energy technique, after integration by parts, we get 

1 



R 



(2.51) 



Taking the real part of this equation and using Poincare 's inequality, we con- 

■-> R 2 — 

elude that H^'H 2 < — t\\F\\ ■ Applying the Poincare's inequality again, we get 

7T 



|^(o,., s )|| 2 <^||i?(o,., s ) 



(2.52) 



which is the desired estimate. 

Using Lemma |2|, we conclude from the three cases above that 



\\{sl-C R )- l f <CR 2 (2.53) 
for \s\ < 2(1 + VI?), C a constant independent of s, R, F, G. 
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3 Conclusions 

The estimates derived for s bounded away from and computed numerically 
for the remaining region of the unstable half plane indicate the L2 norm of the 
resolvent of the operator Cr to be proportional to R. Deriving the estimates 
analytically for the whole unstable half plane is still an open problem as far as 
we know. The method proposed here can be helpful to solve it. It may also 
be possible to adapt this method for the 3 dimensional problem, and with this 
approach it may be possible to give a satisfactory solution to the problem, or 
at least to perform more convincing computations. We hope to address these 
questions in the future. 
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